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We report on the perspective measurements of inclusive particle production in high-energy p-p 
collisions with data to be collected by the LHCb experiment at CERN's LHC. These include VO and 
D meson production studies, which can be based on a minimum bias sample, as well as charmonia 
production studies, which need a muon-triggered samples. Using reconstructed J/V^ ^ decays, 
both the prompt J/'0 and b ^ J/ip production cross-sections will be determined, in the forward 
pseudo-rapidity range of 2-5 covered by LHCb. Due to the large production rate, such analyses will 
be possible with very small integrated luminosities of the order of a few pb~^. Other charmonia 
related measurements will also be discussed, such as that of the J/ip polarization at production or 
of the production of some of the new X, Y and Z states. 



I. INTRODUCTION 

Dedicated to the study of the 6- flavour quark sector, 
the LHCb experiment will take data in proton-proton 
collisions at the CERN's LHC. At an energy in the 
center of mass of 14 TeV, the cross section for bb pair 
production is 500 pb~^, so that, with a nominal lu- 
minosity of 2 • 10^^ cm~^s~^, 10^^ bb pairs will be 
produced in one year (10^ s) of data taking. Within 
this frame the CP violation and rare decays will be 
studied and CKM matrix tests will be performed on 
the full 6-hadrons spectrum as well as in the charm 
sector in search for hints of new physics i2|]. The 
detector of the LHCb experiment is presented in ^ 

It is expected that the 2009 - 2010 LHC data tak- 
ing will start at the center of mass energy = 7 
TeV (may go up to 10 TeV) and an integrated lumi- 
nosity of 100 to 200 pb~l is anticipated. With these 
conditions around 10^ minimum bias events will be 
collected during the first days of data taking. These 
data will allow to perform various production studies 
on VO particles, D mesons and J/V^'s, which will be 
presented in detail in this article. 



II. LHCB DETECTOR 

The LHCb detector is a single-arm spectrometer 
placed in the forward region of the p-p interaction 
point. With an angular coverage, with respect to the 
beam-axis, from 10 (15) to 300 (250) mrad in the 
bending (non-bending) plane, it has an acceptance 
from 1.9 to 4.9 in rapidity, suited to collect the bb 
quarks production which, at the LHC energies, is well 
peaked and correlated in the forward and backward 
region. Here, just a brief description of the LHCb de- 
tector is given, more details can be found in Ref. Q. 

The LHCb detector is divided in several subdetec- 
tors, a simple layout is shown in Fig. [H Around the 
interaction point is the Vertex Locator, a silicon de- 
tector that measures the radial and angular positions 
of charged tracks. The momentum measurement is en- 
sured by the dipole magnet (integrated field of 4 Tm) 
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FIG. 1: Layout of the LHCb detector. 



and the tracking system. The latter is subdivided in a 
Trigger Tracker (TT), a silicon micro-strips detector, 
and three tracking stations (T1-T3) made up of sili- 
con micro-strips for the inner part and straw tubes for 
the outer part. Particle identification, and in particu- 
lar tt — K separation, is ensured by two Ring Imaging 
Cherenckov detectors, RICHl and RICH2. Identifi- 
cation of muons is given by the MUON system, com- 
posed by one detector station (Ml) placed upstream 
of the calorimeter system and 4 downstream (M2-M5); 
the stations are build up from MWPC's with the ex- 
ception of the very inner part of Ml where triple-GEM 
detectors are exploited. Finally, energy measurement 
is made by the calorimeter system: a Scintillator Pad 
Detector (SPD) and Pre-Shower (PS), the shashlik 
Electromagnetic Calorimeter (ECAL) and a hadronic 
calorimeter (HCAL) with Fe and scintillator tiles. A 
summary of the expected performances of LHCb is 
shown in Table H 

A fundamental feature of LHCb is its trigger sys- 
tem. The rate reduction from 40 MHz, LHC bunch 
crossing frequency, to 2 kHz, at which events are writ- 
ten on tape for later analysis, is done by two trigger 
levels. The Level trigger is hardware based, build 
up of custom made electronics, and reduces the rate 
from 40 to 1 MHz mainly requiring particles with high 
transverse momentum {pt) respect to the beam direc- 
tion in the calorimeters and muon system. The High 
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Description 



TABLE I: Resume of the LHCb detector performances. 

Performance 



Momentum resolution 
Energy resolution (ECAL) 
Energy resolution (HCAL) 
5-hadrons mass resolution 
Primary [secondary] vertex position 
Impact parameter 

Time resolution on 6-liadrons proper lifetime 
Kaon identification 
Muon identification 



a(p)/p-0.4% 

a{M) - 14 MeV/c^ 
cr(x) - 50[150] /im 

(j{IP) 14 + 35/p^(GeV/c) /Xm 
cr(t) - 40 fs 

s{K) - 95% at 5% of tt/K mis-id. 
- 94% at 3% of TT, K/iJ. mis-id. 



Level Trigger will exploit the full event data from the 
detector to select events at 2 kHz. It is software based 
and its algorithms will evolve with the knowledge of 
the apparatus performance and physics programme 
leading to high flexibility. 



III. PHYSICS WITH MINIMUM BIAS 

As already said the LHC 2009-2010 run conditions 
will not be the nominal ones, nevertheless, as soon as 
the collisions start, lot of interesting physics will be 
available to be studied. Already within the minimum 
bias events, many measurements can be done: in par- 
ticular in Fig. [2] is shown the physics reach as a func- 
tion of the number of minimum bias collected. The 
physics reach is defined as the relative cross-section 
of a process (with respect to the minimum bias one) 
times the global efficiency for LHCb to detect it. As 
it can be seen and A production can start to be 
studied with about 10^ minimum bias events {i.e. sec- 
onds of data taking); this study is described in ^IVi 
With 10^ minimum bias D production can be stud- 
ied, as described in ^Vl as well as Jip production 
(iVTj), while for b ^ J/t/j something more than 10^ 
minimum bias events will be needed. This last sam- 
ple, 10^ events, is taken as a reference for the per- 
formances of the analysis studies presented in the fol- 
lowing. It is to be bared in mind that such a sample 
corresponds to days of data taking at nominal condi- 
tions, so these studies will be performed as soon as the 
proton collisions will take place; in particular the men- 
tioned physics channels, apart from being interesting 
by themselves, are propaedeutic for the main analyses 
of LHCb, helping either to understand the detector 
performances and as building blocks of 6-hadrons de- 
cays {e.g. B ^ Ji/jKg). 



IV. VO PRODUCTION 

The hadronization process is still not well under- 
stood within the theoretical frame. At hadron col- 
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FIG. 2: Physics reach at LHCb as a function of the number 
of minimum bias events collected. For each channel the 
relative cross-section with respect to minimum bias times 
the efficiency is shown. The straight line corresponds to 
100 events for each channel. 



liders the strange quark production can give hints on 
the hadronization process since strange quarks are not 
present in the initial state valence quarks. Various 
phenomenological models are available, which are of- 
ten different as philosophy, that can reproduce the 
fragmentation into strange hadrons and explain phe- 
nomena and distributions. Within Monte Carlo event 
simulators (e.g. PYTHIA) different tunings are avail- 
able; but, although they do agree at Tevatron energies, 
there are divergences in the extrapolation at the LHC 
energies. Hence the study of VO production at LHC 
(namely A's and Kg) can shed some light in these 
processes. 

There are also experimental motivations that justify 
VO production studies. First of all these decays need 
minimal requirements from the detector point of view. 
Just information from the VELO and the Tracker is 
needed; moreover the simple minimum bias trigger 
will be available from the beginning of data taking. 
Being A's and Kg masses well known, crosschecks of 
the momentum calibration can be made against in- 
variant masses. No particle identification is required 
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so that from these samples RICH cahbration wih be 
available as well. Finally these particles are building 
blocks of subsequent and more complex analysis. 



A. VO analysis 

The study of VO particles in LHCb will start from 
the following channels: K^g tt+tt", A P7r~, 
A p7r~^. The results that will be presented in the 
following come from a study based on simulated min- 
imum bias events with full detector simulation. 

In order to deal with clean events, the ones with just 
one primary vertex have been selected. The analysis 
starts from combining two opposite charged tracks, in 
particular just long tracks are used: i.e. tracks with 
hits along the whole Tracker and in the VELO. As 
already said, no particle identification information is 
required, the analysis being based only on geometry 
and kinematics. 

The Kg analysis requires just the following two 
simple cuts: one on the distance of closest approach 
(DOCA) between the two tracks (required to be less 
than 0.2 mm) and one on the Kg proper time cr > 
4 mm. The A selection imposes the same cuts, plus a 
cut on the impact parameter of the A with respect to 
the Primary Vertex (PV), which is required to be less 
than 0.1 mm. In order to distinguish between Kg, 
A and A decays, without particle identification, the 
Armenteros-Podolanski plot will be used: the trans- 
verse momentum (pt) of decay products with respect 
to the mother particle is plotted versus the longitudi- 
nal momentum asymmetry, which is defined as 



a 



■Pl 



Pt 



-Pl 



where pf^ is the longitudinal momentum (with respect 
to the mother's direction) of the daughters particles. 
The Armenteros-Podolanski plot, for the studied sam- 
ples, is shown in Fig. [3] where the three signal regions 
can be easily distinguished one from each other and 
from the combinatorial background in the lower part 
of the plot. Giving the possibility to separate the three 
mentioned decays, this technique leads to kinematic 
particle identification so that unbiased pions and pro- 
tons samples are available for PID (and hence RICH) 
calibration. 

After this very simple analysis, the invariant mass 
distributions for Ks and A's can show already very 
prominent peaks with good signal to background ra- 
tios; the distributions are shown in Fig. [H The pro- 
duction measurement of VO particles will be made in 
Pt and pseudorapidity (77) bins. Background subtrac- 
tion will be calculated for each bin using the invari- 
ant mass distribution sidebands fit. Within a typical 
sample of 10^ minimum bias events around 2-10^ re- 
constructed Kg are expected and around 10^ A and 
A. 




FIG. 3: Armenteros-Podolanski plot for minimum bias 
events selected as described in ^IV Al Typical Kg , A and A 
ellipses can be seen, apart from some combinatorial back- 
ground in the lower part of the plot. 



Apart from production studies, also the baryonic 
asymmetry will be measurable in the A sector. The 
A/ A ratio, in fact, provides discrimination between 
different hadronization processes [1]. In Figure [5] (a) 
the theoretical expectation for A/ A as a function of 
77 is shown as calculated with different tunings of the 
Monte Carlo generators. As it can be seen from the 
plot, the LHCb experiment will have a larger sensitiv- 
ity with respect to other LHC experiments due to its 
coverage in the forward region of p-p interactions. A 
first analysis within the LHCb experiment, based just 
on 480k Monte Carlo minimum bias events leads to 
the plot shown in Figure [5] (b): the relative error on 
the determination of A/ A ratio is shown as a function 
of the pseudorapidity. As it can be seen the error is al- 
ready at the level of 20-30% which, extrapolated to the 
reference sample of 10^ minimum bias events, leads to 
expected statistical errors at the level of 1.5% which 
will provide discrimination between different models. 



V. D MESON PRODUCTION STUDIES 

During the first days of data taking, as well as 
strange production, charm production will also be 
studied. In particular early studies on D mesons pro- 
duction will be made, by means of the following de- 
cays: K~7r~^ and iiT" 7r+7r+ plus charge 
conjugate modes. The analysis strategy for the se- 
lection of these decays is very similar to the one for 
VO searches. Again only VELO and Tracker infor- 
mation are required, no particle identification being 
needed from RICH, and using only kinematic and ge- 
ometric variables. No use of significance variables {i.e. 
variables divided by their errors) will be made in or- 
der to avoid biases coming from errors not properly 
understood during the early stage of the experiment. 
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FIG. 4: Invariant mass distributions for K% (a) and A (b) 
decays after the selection described in ^IV Al 



Moreover particle ratios will only be studied in order 
to cancel out largely systematics. 
Finally a use of Multivariate analysis techniques is 
made in order to reduce background. The particular 
algorithm used in this case is called Ripper, which is a 
rule based classifier. Analysis studies have been based 
on 9.5 M Monte Carlo minimum bias events. Taking 
the K~^7r~ decay as an example, the following 

variables have been considered both for a traditional 
cut based analysis and for a Multivariate analysis: pr 
of the daughter particles, pr of the D^, impact pa- 
rameters of daughters with respect to PV and angle 
between the two IP vectors, flight length, DOCA, 
impact parameters of D^. The results of the analy- 
sis can be seen in Fig. [6l where the invariant mass 
distributions obtained by selecting events with the 
two techniques are shown: by setting the Multivariate 
method to have the same efficiency on the signal as 
the cut based method, the background is reduced by 
the former of a factor ~ 3. 

From the analyzed sample an average of about 200 
events for each species has been, for the expected anal- 
ysis sample of 10^ minimum bias events around 2000 
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FIG. 5: Barionic asymmetry: (a) A/ A ratio as obtained 
from different tunings of the Monte Carlo generator; it can 
be seen that the LHCb coverage is in a more sensitive re- 
gion in respect to ATLAS and CMS; (b) Plot of the relative 
error on the A/ A ratio, as a function of the pseudorapidity, 
as obtained using 480k simulated minimum bias events. 



events for each species are foreseen. This results lead 
to an expected error on the particle ratios in the sen- 
sitivity region of cr(^) = 5% and cr(^) = 6%. 

VI. J/V^ PRODUCTION 

The production of prompt J/ip is not com- 
pletely understood. Using the Non-Relativistic QCD 
(NRQCD) with Colour Octet Model is possible to 
reproduce the transverse momentum spectrum mea- 
sured at Tevatron; unfortunately this model foresees 
also an increasing of the transverse polarization with 
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(a) DO invariant mass: cut-based analysis 



A. J/ijj selection 

The easiest J/^Ij decay to look for in order to mea- 
sure its production is the one into two muons. Even 
if, as already said, many J/^jj will be already present 
in the minimum bias sample, the use of the Muon 
trigger, which is one of the main blocks of the LHCb 
trigger, will allow to have many more events. The 
selection of J/ijj candidates starts from two opposite 
long tracks with hits also in the Muon stations. A cut 
is made on the of the track and of the vertex of 
the two tracks. One of the two muons is required to 
have a transverse momentum greater than 1.5 GeV. 
Moreover a cut on the identification likelihood (L) is 
applied, with an efficiency on the muons ^ 90% at 
a level of pion mis-identification of ^ 1.4%. 

About 8% of the J/^l) produced at LHC are expected 
to come from 6-hadron decays; in order to distinguish 
between these and the prompt J/ip^s the following 
variable will be used: 
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(b)DO invariant mass: MVA method 

FIG. 6: DO invariant masses distributions as obtained from 
a standard cut-based analysis (a) and from a Multivariate 
analysis (b). In both, in red is shown the background dis- 
tribution and in black the sum of signal and background. 



, _ ^J/iIj — ZPV 

^ - JhP ^^/^ 

Pz 

where zpy and zj/^ are the two vertices coordinates 

and pi^'^ the momentum both along z and mj/^ the 
mass of the J/V^. From studies at the generator level 
t is known to be a good approximation of the 6-quark 
lifetime. In particular in Fig. [71 the distribution of t 
is shown together with the 6-decay time and the two 
distributions are close one to each other. In Fig. [8] 
instead, the distribution of the t variable for different 
J/^Ij contributions is shown. In particular, apart from 
a ffat combinatorial background (the shape of which 
will be estimated from J/ip invariant mass sidebands), 
the prompt J/^l) peak at t = and the J/ijj from 
h exponential tail can be distinguished. Within the 
long tail of the distribution there is also a component 
due to true prompt J/?/^ associated to the wrong pri- 
mary vertex giving so a wrong t value. The amount 
of this component will be estimated by using true re- 
constructed J/V^ and associating them with primary 
vertices coming from other events chosen randomly. 



B. J/V^ measurement 



the pt-, which has not been seen at Tevatron. Alter- 
native models {e.g. Colour Evaporation) also predict 
the same pt spectrum but cannot account for the po- 
larization. Hence cross-section and polarization mea- 
surements are important to understand charmonium 
production. Given the big cross-section of J/^l) pro- 
duction and the unique LHCb coverage, large data 
samples will be available already in the first days of 
LHC running._ At the same time also the production 
of J/ij) from hh decays will be measured. 



Using the first 5 pb~^ of integrated luminosity 3 
million reconstructed J/?/^ are expected. The event 
counting will be made in pseudorapidity and trans- 
verse momentum bins. In order to compute the num- 
ber of J/il) in each bin a fit of the invariant mass dis- 
tribution will be made. In particular the plot shown in 
Fig.[9]shows this distribution obtained by just using 19 
millions minimum bias events (corresponding to sec- 
onds of data taking at nominal luminosity), the mass 
resolution obtained is about 11 MeV. Monte Carlo 
data will be used to correct for detector acceptance 
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FIG. 7: h quark lifetime distribution studied at the event 
generator level in blue and t variable distribution in red; 
it can be seen that the two distribution are very similar. 
Prompt Jijj events in this plot would lie on the very far 
left. 
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FIG. 9: J/V^ invariant mass distribution (in MeV) as ob- 
tained by applying the selection sketched in WI Al to a 
sample of 19 millions minimum bias events. 
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FIG. 8: Distribution of the t variable for different com- 
ponents: the blue curve is the prompt Jpsi distribution, 
the red curve are Jip from b, and the magenta the sum of 
the two. Within the blue curve, the long tail is due to the 
combination of true Jip with wrong Primary Vertices. 

and efficiencies for trigger and offiine selection. Fi- 
nally in order to measure the absolute cross-section, 
integrated luminosity measurements will be needed. 

As already said, an important parameter to mea- 
sure will be the J/ip polarization 0; this is defined as 
the angle between the /i+ momentum in the J/tJj rest 
frame and the J/t/j momentum in the lab frame. In 
the helicity frame the distribution is parametrized as 

-— cx 1 + a COS [0) 

acos{0) 

where a = 1 means transverse polarization, a = —1 
longitudinal polarization and a = no polarization. 
Simulation studies have shown that the LHCb geome- 
try and LO Trigger efficiency will induce a fake polar- 
ization on the J/ip reconstructed sample, even if the 
original sample were not polarized; so that in order to 
measure the J/i/j cross-section, its polarization should 
be taken into account. Polarization measurement then 



will be made and will be also used to discriminate be- 
tween different production models. 



VII. Xc PRODUCTION 

From Tevatron measurements it's known that about 
30% of J/i/j come from Xc(i,2) ^ J/i^l and hence 
have different polarization. The production of Xc(i,2) 
is also interesting by itself being the ratio R^^ = 
o'{Xi)/^{X2) important to distinguish different pro- 
duction models. While the NRQCD predicts R^^ to 
be close to 1, Colour Evaporation Model predicts it to 
be 5/3 from the expectation from spin-counting. 

In order to select Xc candidates a photon with 
Pt > 500 MeV is combined with a J/?/; selected as 
before (restricted to a ±40 MeV mass window). The 
distribution of AM = M^^^ — M^^ is then used to ex- 
tract the Xci and Xc2 signals. This distribution, shown 
in Fig. \W\ will be fitted with two gaussians which 
represent the signal and a background component 
parametrized as P(AM) = AM^o •exp(-ci • AM-C2 • 
AM^). The obtained resolution is a{AM) ^ 27 MeV 
which has to be compared with the mass difference of 
the two resonances, m(xc2) — ^(Xci) = 46 MeV; the 
achievable separation is limited but still some sensi- 
tivity is present leading to a discrimination between 
Xci and Xc2 and to the measurement of R^^ . 



VIII. OTHER QUARKONIA 
MEASUREMENTS 

Beyond the already mentioned studies, measure- 
ments will be performed also on other charmonia and 
bottomonia channels. The V^(26') production and in 
particular the o-{ip{2S))/o-{J/tlj) ratio will be stud- 
ied. Moreover, similar measurements to the char- 
monia case will be performed also for bottomonium 
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A M/GeV 

FIG. 10: Distribution of the AM, as defined in ^Vlll for 
inclusive J/'0 events. The AM for Xc(i,2) ~^ J 1^1 events 
peak can be seen. The red and green curves are the gaus- 
sian distribution used to represent the signal; the blue 
curve represents the background (J/'0 not from Xc) dis- 
tribution. 

production and spectroscopy {e.g. T{1S^2S^3P) res- 
onances, Xb ^l)' 

The exotic X, Y and Z charmonia states, recently 



observed, will also be studied at LHCb. In particu- 
lar is under investigation the possibility to measure 
the X(3872) quantum numbers by means of an angu- 
lar analysis of the decay X(3872) J/?/^7r+7r~; this 
analysis will be even more sensitive if studied in the 
frame of the XK^ decay where the X is ex- 

pected to be produced polarized: in this context will 
be possible to distinguish if the X has J^^ equal to 
1++ or 2-+. 



IX. CONCLUSIONS 

To summarize, the LHCb experiment, ready to take 
data at the LHC proton collisions, will be able to ex- 
ploit the very first data, apart for detector calibrations 
and tunings, also for very interesting physics analy- 
ses. Within a sample of 10^ minimum bias events it 
will be possible to probe hadronization models with 
VO production studies and D-mesons ratios will be 
measurable with about ~ 5% error. Quarkonia spec- 
troscopy will be also investigated and J/ i/j production 
and polarization will be measured. 
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